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Abstract 

We comprehensively analyse the theoretical prediction for the Lamb shift in muonic 
hydrogen, and the associated determination of the proton radius. We use effective 
field theories. This allows us to relate the proton radius with well-defined objects in 
quantum field theory, eliminating unnecessary model dependence. The use of effective 
field theories also helps us to organize the computation so that we can clearly state 
the parametric accuracy of the result. In this paper we review all (and check several 
of) the contributions to the energy shift of order a 5 , as well as those that scales like 
a 6 x logarithms in the context of non-relativistic effective field theories of QED. 
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1 Introduction 

The measurement DIZI of the Lamb shift in muonic hydrogen, 


E(2P 3/2 ) - E(2S 1/2 ) = A E e * p = 202.3706(23) meV (1.1) 


and the associated determination of the root mean square electric radius of the proton: 
r p = ^ j(r p ) = 0.84087(39) fin has led to a lot of controversy. The reason is that this number 
is 7.1a away from the CODATA value, r p = 0.8775(51) fm [3j. This last number is an average 
of determinations coming from hydrogen spectroscopy and electron-proton scattering^. In 
order to asses the significance of the discrepancy, it is of fundamental importance to perform 
the computation of the Lamb shift in muonic hydrogen (in particular of the errors) in a 
model independent way. This was done in Ref. [6]. In that Letter we revisited the theoretical 
derivation of the Lamb shift using effective field theories (EFTs) and obtained the following 
expression 


A E? 


206.070(13) — 5.2270(7)—A 

mb 


meV. 


( 1 . 2 ) 


Using this result and Eq. (jl.lj) we then obtained 


r p = 0.8413(15) fm, 


(1.3) 


which is at 6.8a variance with respect to the CODATA value. Therefore, the proton radius 
puzzle survived our model independent analysis. The good point now is that the EFT 
analysis allows us to have a parametric control of the uncertainties, which are of the order 

3 

of uncomputed terms of m p a e ). This parametric control of the uncertainties 

allowed us to obtain a model independent estimate of the error, which is dominated by 
hadronic effects. 

EFTs help organizing the computation by providing with power counting rules that asses 
the importance of the different contributions. This is specially so for the muonic hydrogen, 
as its dynamics is characterized by several scales: 

m p ~ m„, rrin ~ m, ~ m r = -, m r a ~ m e . 

m p + m p 

By considering ratios between them, the main expansion parameters are obtained: 

m-rr m u 1 m e m r a m r a 1 2 1 

— ~ —- ta — , — ~-~-~ a ~-. (1.4) 

m p m p 9 m r m r m r a 137 

For our evaluation we used potential non-relativistic QED (pNRQED) [7]. Particularly 
relevant for us is Ref. [8], which contains detailed information on the application of pNRQED 
to the muonic hydrogen. Since pNRQED describes degrees of freedom with E ~ m^a 2 , any 
other degree of freedom with larger energy is integrated out. This implies treating the proton 

1 The latter though has been challenged in Refs. 110 , and its exclusion would certainly diminish this 

tension. 
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and muon in a non-relativistic fashion and integrating out pions (and Delta particles). This is 
the step of going from ffeavy Baryon Effective Theory (HBET) [9] to Non-Relativistic QED 
(NRQED) [10], By integrating out the scale m p a, pNRQED is obtained and the potentials 
appear. Schematically the path followed is the following (A = tua — m p ): 

HBChPT ^ NRQED W pNRQED. 

A detailed explanation of the matching computation between ffBET and NRQED was given 
in Ref. [TTj . This corresponds to the hadronic part of the computation presented in Ref. [6]. 
It is one of the main motivations of this paper to give the details of QED-related part of the 
analysis in Ref. [6) . This means to analyse the potentials that contribute to the given order, 
as well as to actually compute the associated energy shifts associated to the potentials and 
the nltrasoft photons. We have made some effort to present the result assuming an arbitrary 
charge for the muon and proton, so that the results can be of use in a more general situation, 
in particular for muonic atoms. This is so because the expressions of the potentials would be 
equal for light muonic atoms after appropriately changing the NRQED Wilson coefficients 
produced by the hadronic scales. Therefore, we will present some results in terms of Z p {= 1), 
Z p {= 1) and Z = Z p Z p (= 1). We also expect that the analysis presented in this paper will 
set the basis for higher order computations using EFTs. 


2 NRQED (fip) 


In the muon-proton sector, by integrating out the ~ m p scale, an EFT for non-relativistic 
muons and protons, relativistic electrons and photons appears. In principle, we should also 
consider neutrons but they play no role at the precision we aim. The effective theory has a 
hard cut-off v <C ro, and therefore pion and Delta particles have been integrated out. The 
effective Lagrangian reads 


A 


= + C e + £( nr) + C N + C 


C 


(NR) 


-NRQED(/P ■‘-'7 t >—e t p ~r ~r ^Ne ~r 

The pure photon sector is approximated by the following Lagrangian 


£ 7 = 


. _ TP» V TP 
4 


d 


O) 


m f 


+ 


do 


mi 



F pv D l F^ v , 


( 2 , 1 ) 


( 2 , 2 ) 


and d^ are generated by the vacuum polarization loops with only muons and taus 
respectively. At 0(a) they read 

d “ = flr + 0( “ 2) ’ 4T, = 6S£ 0( “ 2) ' (2 ' 3) 

The hadronic effects of the vacuum polarization are encoded in d 2 ; 

<k = ^n;,(0) = + df* + 0(a 2 ). (2.4) 
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n;(0) is the derivative of the hadronic vacuum polarization (we have defined 11^ (—k 2 ) = 
—k 2 n/(0) + ...). The experimental figure for the total hadronic contribution reads 11^ ~ 
9.3 x 10 -3 GeV~ 2 [12]. Following standard practice, we have singled out the contribution due 
to the loops of protons (assuming them to be point-like) in the second equality of Eq. (12.4p . 
Note though that is still of order a. 

The electron sector reads ( iD At = id M — eA 

C e = l e {ilp -mi e )l e . (2.5) 

We do not include the term 

- , (2.6) 
mf, 

since the coefficient gi e is suppressed by powers of a and the mass of the lepton. Therefore, it 
would give contributions beyond the accuracy we aim. In any case, any eventual contribution 
would be absorbed in a low energy constant. 

The muonic sector reads 


£(NR) 


i&i + 


2m,, 


+ 


8 ml 

h 1 


+ e 


r (/d 

l f 

2 m,, 


a B 


+e 


Jp) 

c p 

8m 2 

A * 1 


[V-El 


le 


r 0d 

c s 

8 ml 

A 6 


(T ■ (D,, xB-ExD, 



(2.7) 


with the following dehnitions: 


iD ° = id n 


- Z.eA 0 , iD, 


Wilson coefficients can be computed order by order in a. 
the fact that dl 3 = 2— Z, 


- zV + Z^e A and Z M = 1. The 
They read (where we have used 


cf = (l + ?£ + 0(a 2 )) , 

cf = (1 + ^ + °(“ 2 )) • 


( 2 . 8 ) 

(2.9) 


Taking the values of the form factors for the muon-electron difference computed in na 
and those for the electron computed in [15], we can deduce the following expression for the 
c o Ms( ,y ) Wilson coefficient!!: 

2 In NRQEDffip), the electron has not been integrated out. Therefore, Eq. (12.101) is not the Wilson 

coefficient of NRQED(^rp). Eq. (12.101) will show up after lowering the muon energy cut-off below the electron 

mass in pNRQED. Still we choose to present it here as otherwise we would be forced to do an extra inter¬ 

mediate matching computation that it is unnecessary to obtain the final result. Since we have integrated 

out the electron, note also that a = 1/137.14... in this equation, i.e. any running associated to the electron 
is written explicitly in Eq. (12.1011 . 
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( 2 . 10 ) 


„(!•)_ 

D,MS 


[V) 


4a n, 

1 H- Z~ In 

3vr ** 


m. 


v* 




7T 


-InM^ 

m P 


+zl 


T( 181n(2) -f 


O (a 3 ) . 


40 , (m u \ 85 4n 2 
27 UJ 81 27 


1^23 _ 9 ' 

324 2^ } 


+ o 




m h 


Note that written in this way one can easily read the 0(a 2 ) Cj and CfT F ni (for the case 
of massive rp quarks) coefficients that would appear in the analogous Wilson coefficient c F 
in QCD. The second line in Eq. (I2.10p would correspond to the CfTpUi term and the third 
line to the Cj one. 

For the Lamb shift computation we perform in this paper we only need c P with 0(a 2 x In) 
accuracy. We also include the finite piece for completeness but neglect 0[m e /m p ) terms. 
Note that analogous 0(a 2 ) terms (changing by m p and either keeping m e or changing it 
by m M ) would exist for if computing the Wilson coefficient as if the proton were point¬ 
like at the m p scale. Even if these effects are small, they should be taken into account for 
eventual comparisons with lattice where typically only the hadronic correction is computed. 

For the proton sector we have 


Cn = 



— e- 


Sp) 

'F 

2 m r . 


a B 


8mj 


cr ■ (D p x E — E x Dp) >Nj 


p > 


( 2 , 11 ) 


where iDjl = id 0 + Z p eA°, iD p = zV — Z p e A and for the proton Z p = 1. The proton Wilson 
coefficients are hadronic, non perturbative quantities. In some cases they can be directly 
related with low energy constants, for instance with the anomalous magnetic moment of the 
proton, k p = 1.792847356(23) [IS]: 

Z 3 a 

= Z p + k p = Z p + Kp ad -| dh—b 0(a 2 ), (2-12) 

Z 3 a 

c { p = Z p + 2n p = Z p + 2^ ad + -5- + 0(a 2 ). (2.13) 

7r 


Note that k p includes O(a) effects. In principle, this is also so for Kp ad , to which we have 
subtracted the proton-associated point-like contribution to the anomalous magnetic moment 
(note that the point-like result is a bad approximation for cp , even though it gives the right 
order of magnitude). The case of cp is more complicated (a more detailed discussion can 
be found in Ref. 0 )- It can be written in the following way in terms of the electromagnetic 
current form factors at zero momentum (iq(0) = Z p ): 

cp ] = Z P + 2F 2 (0) + 8F((0) = Z P + 8 m 2 p dG Pf[ q ^ . (2.14) 

d( l q2=0 
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This object is infrared divergent, which makes it scale and scheme dependent. This is not a 
problem from the EFT point of view but makes the definition of the proton radius ambiguous. 
The standard practice is to make explicit the proton-associated point-like contributions to 
the computation. In practice this means that one uses the following definition for the proton 
radius 


r (p) _ 

D,MS 


v 


= Zp +-3^ In 
P 3 7T 




r p m p 


+ 0(a 2 4 


(2.15) 


In other words (up to 0(a 2 ) corrections) 


(p)_ 
D, MS 


\m. 


v — 2 2 

Z p = —r p m p . 


(2.16) 


Note that r p includes 0(a) terms in its definition. This should be kept in mind when 
comparing with lattice determinations. Note also that it is not natural to set v = m p , or, 
in other words, to assume that the proton is point-like up to (and beyond) the scales of 
the proton mass; | r 2 m 2 ~ 21.3, to be compared with ”1” for a point-like particle. This 
illustrates that the point-like result does not even give the right order of magnitude of 
£jv e refers to the four-fermion operator made of nucleons and (relativistic) electrons. It 
will not contribute to the spectrum at 0(m r a 5 ). Therefore, we will not consider it any 
further. For a more detailed discussion see Ref. [TTl . 

Finally, we consider the four-fermion operators 


% - ZT2 N > N V • ( 2 - 17 ) 

" L p " L p 

Again in this case it is common practice to single-out the proton-associated point-like contri¬ 
bution. Note that this assumes that one can treat the proton as point-like at energies of the 
order of the proton mass. We have already seen that this is a bad approximation for cd and 
other Wilson coefficients. Nevertheless, we keep this procedure for the sake of comparison. 
Therefore, 


TY1 

C 3 (^) =- -d s {v) + C 3 ad + 0(a 3 ), 

(2.18) 


c 4 = ~^d v + cf 1 + 0(a 3 ), 

(2.19) 


m p 


where the point-like Wilson coefficients read as follows: 


d s (v) 


(ly 


Z 2 a 2 


m 2 — m 2 


p l 


m, 


Z 2 a 2 


m f 


m. 


2 m /1 m p 


In —y 


m 

In —h. 

m 2 


+ 



( 2 . 20 ) 

( 2 . 21 ) 


3 Although not directly relevant for the specific computation of this paper, note that this also happens 
for the Wilson coefficients and ca 2 (for the definition see Ref. [ 8 ]), for which their physical values are far 
from zero: ca 2 — 12 and ca 2 — —72, even though for a point-like particle their values would be ”1” and ”0” 
respectively (up to 0(a ) corrections). 

4 The coefficients C 3 and C 4 should actually read c 3 ^ and , as they actually depend on the nucleon and 

lepton the four-fermion operator is made of. Nevertheless, to ease the notation we eliminate those indices. 
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The expression of d s should be understood in the MS scheme, d v on the other hand is finite. 
d s was computed in Ref. |l8j and d v in Ref. [TO] . 

Cg ad encodes all the hadronic effects to the spin-independent four-fermion Wilson coeffi¬ 
cient. At 0(a 2 ) it is generated by the two-photon exchange contribution. Since C 3 ad depends 
linearly on the muon mass, it is dominated by the infrared dynamics and diverges linearly 
in the chiral limit. This produces an extra m^/m^ suppression with respect to its natural 
size, and allows us to compute the leading pure-chiral and Delta-related effects in a model 
independent way. The complete matching computation between HBET and NRQED was 
made in Ref. [lTj to which we refer for details (partial results can be found in [8, T9] , and in 
Ref. [2D] in the context of relativistic baryon effective theory). Overall we obtained 


had 

c 3 


a 


m. 


1 + #X + '" 


+ 0 a 


2 "V 


A, 


QCD 


= a 


2 m H 


47.2(23.6) (tt), 

56.7(20.6) (tt + A), 


( 2 . 22 ) 


where the upper and lower numbers refer to the matching computation with only pions, or 
with pions and the Delta particle, respectively. For comparison, the value c| ad = a! 2 ^54.4(3.3), 
which follows from the analysis in Ref. [2T], was used in Ref. [2]. We refer to Ref. mi for a 
detailed discussion on the status of these hadronic determinations and focus on the QED-likc 
computations in this paper. 

c|j ad encodes all the hadronic effects to the spin-dependent four-fermion Wilson coeffi¬ 
cients. As in the previous case, this coefficient diverges in the chiral limit. Nevertheless, it 
only does so logarithmically (unlike in the previous case, where the divergence was linear). 
Such computation can be found in Ref. mi. Still it is possible to determine from the 
analogous one of the proton-electron four-fermion operator determined in Ref. [17] . This 
was done in Ref. m from where we quote the result 

cf d _ _ 46a 2 (2.23) 


3 pNRQED 


After integrating out scales of O^m^a ~ m e ), the resulting effective theory is pNRQED [7]. 
This EFT naturally gives a Schrodinger-like formulation of the bound-state problem but still 
keeping the quantum field theory nature of the interaction with ultrasoft photons, as well 
as keeping the information due to high energy modes (of a quantum field theory nature) in 
the Wilson coefficients of the theory. pNRQED has been applied to hydrogen [22], positron- 
ium [ 23 ] and muonic hydrogen mm providing with much of the information needed for this 
paper. In particular in the last reference the explicit form of the Lagrangian was presented 
(up to 0(m r a 5 )). We repeat it here but generalized to the case of arbitrary charges: 


^pNRQED 


= /#xd3x S t(x.X. i ){ iao -|7 + £ 


+ 


p 


8 mf, 8 2 M 


(3.1) 


-V(x,p, <Ti,CT 2 ) 


Z p m p + Z v m 


p' 


m p + m p 


x • E(X, t) ^(x, X, t) — I d 3 X-F pu F^, 
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where M = m u + m„, m r = m,L ™‘ v , x and X, and p and P are the relative and center of 
mass coordinate and momentum respectively. 

V can be written as an expansion in 1/m At , l/m p , a, ... We will assume 1/r ~ m e (which 
is realistic for the case at hand) and that m M <C m p . We then organize the potential as an 
expansion in l/m p : 

T(x, p, (Ti, (T 2 ) = H (0) (r) + H (1) (r) + P (2) (r) H - , (3.2) 

where 

P (n) oc —. (3.3) 

ml 

We will also make the expansion in powers of a explicit. This means that 

K J_ a r (3.4) 

ml 

p(O’i) = — — has to be included exactly in the leading order Hamiltonian to yield the 
leading-order solution to the bound-state problem: 


p 2 Za 
2 m r r 


(3.5) 


Thus, the contribution to the energy of a given potential is 


(V M ) ~ m p a 1+n+r 


up to large logarithms or potential suppression factors due to powers of l/m p . Iterations of 
the potential are dealt with using standard quantum mechanics perturbation theory produc¬ 
ing corrections such as: 


{V M ■ ■ ■ V {m ’ s) ) ~ m At a 1+n+r+(1+m+s_2) (3.6) 

and alike. Therefore, in order to reach the desired 0(ma 5 ) accuracy, has to be computed 
up to 0(a 4 ), V^ up to 0(a 3 ), up to 0(a 2 ) and up to O(a). 

3.1 The static potential: 

The Fourier transform of V ^ reads 

V (0) = -4t tZ^^- = ^ ( °’ n) > (3.7) 

n =1 

which in fact defines ccp, the coupling constant associated to the static potential, which is 
gauge invariant. The contribution associated to the electron vacuum polarization (n(0) = 0) 


n {k 2 ) = alliik 2 ) + a 2 U 2 (k 2 ) + a 3 U 3 (k 2 ) + ... 





provides with another very popular definition for the effective coupling that enjoys the nice 
properties of gauge invariance and scheme/scale independence: 

1 2 3 4 

a eS (k) = a TTTu~^ = «- —-Hi + ^(n? - n 2 ) + ^3 (-n?+2n 1 n 2 - n 3 )+ o(a 5 ). (3.8) 

a e s corresponds to Dyson summation. If we express aty(k) in terms of a eff (/c), we have 
ay(k) = a eS {k) + J2 Z;Z^a^ m) (k) = a eS (k) + 5a(k), 6a(k) = 0(a 4 ). (3.9) 

n,m =0 

n+m=even>0 


D (0,1) = 

k z 


(3.10) 


is nothing but the Coulomb potential. In order to achieve 0(m r a 5 ) accuracy we need to 
know nW, II (2 \ n® and the leading, non-vanishing, contributions to a^°\k), a^ 2 \k) and 


a 


(i,i) 

eff 


(k). 



Figure 1: One-loop electron vacuum polarization contribution to the static potential. 
The next-to-leading order term of the static potential is displayed in fig. 1 and reads 

,ck 2 Ilif—k 2 ) 


Vi°/\k) = 4 tt Z- 


where 


and 


n^jfc 2 ) = k 2 


dq 2 


7T k 2 
1 


q 2 (mlq 2 — k 2 ) 


u(q 2 ), 


= 1 + 1 • 


q 

Thus, we may write for the potential in position space 


K 


( 0 , 2 ), 


Za a 


poo 


dq 2 


VP 


r) =-/ — u(q )e 


2\ -2m e rq 


r IT J 4 q* 


(3.11) 

(3.12) 

(3.13) 

(3.14) 
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Figure 2: Diagrams contributing to V^ 0,3 ). 


The next-to-next-to-leading order term of the static potential is produced by the diagrams 
depicted in Fig. [2j which can be understood as a correction to the vacuum polarization. It 
was computed by Kallen and Sabry m and reads 


Vi°/\k) = 4vrZ 


U' 2 1 (k 2 )-U 2 (k 2 ) = k 2 


a 3 n 2 (-k 2 ) -n 2 (-k 2 


7T Z 


k 2 


dq 2 


q 2 (m 2 q 2 — k 2 ) 


«<V), 


(3.15) 

(3.16) 


where 
u^ 2 \q 2 ) = 


19 55 2 1 4 3 — t 2 , 

t (-+ —r 2 - -r 4 -n 

1 24 72 3 2 


64r 4 


(l-r 2 r 


In 


1 + T 


1 — r 

+ (3 + 2r 2 - r 4 ) ( 2Li 


— H—-r 2 --r 4 -|—r 6 + ( - + t 2 -) In 

16 8 16 6 \ 2 2 1 


(1 + t) ; 

8 r 2 




1 — r 
1 + r 


L4; 


-1 + T 

1 + r 


with 


Li 2 (» = - 


du —-- zeC\[l,oo) and r = \ jl - 


'0 


u 


T 


(3.17) 


(3.18) 


The next-to-next-to-next-to-leading order term of the static potential coming from the 
vacuum polarization reads 


V^ A \k) = 4vrZ 


a 4 -n 3 (-k 2 ) + 2II 1 (-k 2 )n 2 (-k 2 ) - n 3 (-k 2 

7T 3 k 2 


(3.19) 


This object (more specifically n 3 ) has been computed in Ref. [25], see also [26] where the 
complete set of diagrams can be found. 

The remaining next-to-next-to-leading order contribution to the static potential is gen¬ 
erated by diagrams that cannot be completely associated to the vacuum polarization, and 
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is encoded in 5a(k). Its sum is constrained to fulfil n + m = even because of the Furry 
theorem. Each a^ m \k) is also gauge invariant. The leading, non-vanishing, contributions 
are a^f\k) and cx^’ 1 ^). They have an expansion in cx themselves. Since each of 

them is (9(cx 4 ), we can approximate them by its leading order expression, which is produced 
by the light-by-light diagrams displayed in Fig. El This object could be deduced from the 
computation in Ref. [27] (we truncate the a eff expressions to its leading order) 

vS,t\k) = ~ {(Zl + 2>g o) (J0 + Z,Z p a^\k)) , (3.20) 

where we have already used that a^°\k) = a^ 2 \k). 


-I-»- fj, 

i 




v -»- 1 ->- p 


(1 : 1 ) 


(2 : 0 ) 


Figure 3: Light-by-light contribution to the static potential. The first and third diagram are 
the contributions to cx^’ 0 ^) and a^ 2 \k) respectively. The second diagram contributes to 


cx 


(i,i) 

eff 


(*)■ 


3.2 The potential beyond the static limit 

In the matching scheme used in this paper (off-shell in the Coulomb gauge) the 1/m potential 
is zero in QED without light fermions (see Ref. [23]). This is no longer so after the inclusion 
of light fermions (the electron) into the computation. Yet, after inspection of the diagrams 
that may contribute, they would produce, at most, 0(m r a 6 ) corrections to the energy, so 
they will be neglected in the following. 

At order 1/m 2 the expression of the potential in momentum space was obtained in Ref. [8]. 
We summarize its different contributions here. If we switch off the electron vacuum polar¬ 
ization effects, the computation would correspond to the muonium case (or positronium 
for the equal mass case). The relevant diagrams in such situation are presented in Figs. Q] 
and [5] (following the classification of Ref. [23] generalized to the non-equal mass case). 
The one-loop diagrams in Fig. [5] produce the potential quoted in Eq. (I3.22|) (in the MS 
scheme). The tree-level diagrams of Figs. [4] produce the potential quoted in Eq. (I3.2ip with 
ex eff = cx. In order to incorporate the electron vacuum polarization effects at one loop we 
replace cx —> cx e fr — cx — cx 2 II 1 (—k 2 ) in Eq. (I3.2ip . This means including the vacuum po¬ 
larization in the 1st, 2nd, 4th, 5th and 6th diagram in Fig. [4j On top of that one has to 
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-Q-*- P 


P 


■g—*- p 


p 



p 





p 




Figure 4: The non-zero relevant diagrams for the matching at tree level in the Coulomb 
gauge. The dashed and zigzag lines represent the A 0 and A fields respectively, while the 
continuous lines represent the fermion and antifermion fields. For the A 0 the circle is the 
vertex proportional to cd, the square to cs (spin dependent) and the black dot to c? 2 , while for 
A the square is the vertex proportional to cf and the other vertex appears from the covariant 
derivative in the kinetic term. The last diagram is proportional to C3 and C4. The symmetric 
diagrams are not displayed. It corresponds to Eq. \3.21\) 

include the contribution coming from Fig. [HI which appears from the Taylor expansion in 
powers of the transfer energy of the vacuum polarization when doing the matching compu¬ 
tation off-shell (for further details see the discussion in Ref. [8]). It produces the potential 
quoted in Eq. (13.23(1 . Overall, the complete expression for the 1/m^ potential at one loop is 
summarized in Eqs. (13.211) . (I3.22p and (13.231) . which we list now 
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Figure 5: The non-zero relevant diagrams for the matching at one loop in the Coulomb gauge. 
The dashed and zigzag lines represent the A 0 and A fields respectively, while the continuous 
lines represent the fermion and antifermion. The interactions for A are the ones which 
appear from the covariant space derivatives in the kinetic term, while for Aq comes from the 
covariant time derivative. The symmetric diagrams are not displayed. They correspond to 

Eq. mm) . 
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(3.21) 

(3.22) 

(3.23) 


where Si = <Xi/2 is the spin of the particle i. We stress again that Eq. (I3.22j) has been 
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Figure 6: Symbolic representation of the leading correction to the static potential due to the 
Taylor expansion of the electron vacuum polarization in the Coulomb gauge in powers of 
k° = Ei — E[. It corresponds to Eq. \3.23\) . 

obtained in the MS scheme. The sum of these three potentials includes all terms of 0(l/( 2,1 )) 
and £>(l/( 2 ’ 2 )): 

V m = nS+vp + vSlp + Uff-Lii + 0(C ,2 ' 3) ), (3,24) 

3.3 The potential in position space 

The matrix elements of the potentials that appear in the energy shifts are more efficiently 
computed in position space. Therefore, we also write the potentials in position space. In 
this case it is convenient to split the potential in a slightly different way than in momentum 
space. In particular, the vacuum polarization contributions are dealt with in an isolated way. 
We follow the notation of Ref. |8]. The contributions coming from tree-level diagrams read 
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where Sp(r) = —4(Si ■ Sj) + 12(Si • r)(Sj • ?). 
The Fourier transform of Eq. (I3.22|) reads 
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Finally, the contributions associated to the one-loop vacuum polarization reac@ 
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where A = m e q. Therefore, with the precision we aim at, we obtain 

V m = viS+y^ip+try p ) , 1 _ loop +v' 0 H M1 )+o(i/' 2 ’ 3 >) = ^i+vM + o(rwi), (3.29) 

where in the second equality we organize the potential terms according to their powers in 
a. This requires expanding the NRQCD Wilson coefficients in powers of a. The leading 
non-vanishing contribution read^l 
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5 Note that the fourth line can be written in a way that makes the angular momentum structure more 
explicit: 
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(3.27) 
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Ar 


Nevertheless, one has to be careful when dealing with the right-hand-side of the equality, as the first and 
last term are separately divergent (but not their sum). 

6 Strictly speaking there could still be some 0(a) included in Kj) ad with the definition we are using, 
similarly as it happens with the proton radius. 
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For the organization of the computation it is also convenient to split l/*- 2,2 -* in the following 
way: 

^ (2 ’ 2) = K ( o-vp + ^vp 2) - (3.31) 

The first term is the potential if we switch off the interaction with the electrons. This is 
a well defined limit, as it corresponds to the case of muonium. The second term is the 
correction to the potential associated to the one-loop electron vacuum polarization. 

Finally, the 1/m 3 potential, which we directly consider in position space, just comes from 
the Taylor expansion of the relativistic expression of the dispersion relation: 


y(3,o) 


There are no 0(a/m 3 ) terms. 



(3.32) 


4 Muonic hydrogen Lamb shift: E(2P 3 / 2 ) ~ ^(25i/ 2 ) 


In this section we review all (and check several of) the contributions to the energy shift of 
order a 5 , as well as those that scales like a 6 xlogarithms in the context of pNRQED. The 
muonic hydrogen Lamb shift is obtained by the combined use of non-relativistic quantum 
mechanics perturbation theory and perturbative quantum field theory computations (when 
ultrasoft photons show up). As we have definite counting rules to asses the relative impor¬ 
tance of the different terms we know when we can stop computing. The application of this 
program to the muonic hydrogen produces the contributions we use in our analysis, listed 
in Table [H Most of the results were already available in the literature, we have checked 
many. We now discuss them focusing on the novel aspects. Note that, even though most of 
the contributions can be associated to a pure QED calculation, the hadronic effects are also 
included in this computation. Their effects are included in the NRQED Wilson coefficients 
discussed in Sec. EJ and are encoded in the different terms of the potential in the Lagrangian 
of pNRQED discussed in Sec. El 

In order to carry out the computations of this paper we use the most updated PDG 
values [T6j for the masses and fine structure constant 

m e = 0.510998928(11) MeV, 

= 105.6583715(35) MeV, 
a = 1/137.035999074(44), 
m p = 938.272046(21) MeV, 

m T = 1776.82(16) MeV. (4.1) 


These numbers update the values used in Ref. [6]. The effect is very small but it changes 
the last digit of our numbers in some cases after rounding. This happens in the second term 
in Eqs. (II. 2115. ip . in Eq. (1 1.3 j) . and in the i) and vii) entries in the table [0 

p n i is the non-relativistic charge density of the nl state. For the n = 2 Lamb shift we will 
need their difference between the S- and P- wave bound state: 


P = P 2 P ~ P 2 S = ( m r Zaf ,2 e 


—m r Zar 


12 


(m r Zar) 2 — ( 1 — 


m r Zar 


(4.2) 
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We will use the following notation: 


SE^ = (nlj\V\nlj) 


and 

SEY = (2P 1/2 \V\2P 1/2 ) - (2S 1/2 \V\2S 1/2 ) 

for the correction to the Lamb shift of a generic potential V. 


V V 


(4.3) 

(4.4) 


Figure 7: 2nd order perturbation theory of the bound-state Green function generated by a 
generic potential V. 


We will represent the 2nd order perturbation theory correction to the bound-state Green 
function generated by a generic potential V by Fig [TJ where the double line represents the 
bound state and the vertices (local in time) the potentials. In case we want to obtain the 
associated energy shift we will compute objects like (and analogous expressions in case of 
different potentials (including permutations)) 


8E% V = (M v by VM 

= j dr 2 dr 1 V>i i (r 2 )V'(r 2 )G^(r 1 ,r 2 )V'(ri)i/> nW (r 1 ), (4.5) 

where 

{E n - h)' = ~ E — En) 1 ^ 4 ' 6 ^ 

G' nl (r 1 ,r 2 ) = (ri| 1 |r 2 ) = lim far x , r 2 ; £) - "j , (4.7) 

(■ Eni-hy E=,E n \ E — E n J 

-0nz( r ) is the bound state wave function of the (nZ)-state and E n is the energy of the state, 
and G(r-i,r 2 \ E) is the Coulomb Green function. 

In order to perform the computation it is specially useful to use the following represen¬ 
tation for a negative energy E = — mr .f x2 a of the Coulomb Green function (see for instance, 
the appendix of Ref. [28]): 


G(r 1 ,r 2l E) = 


mlZa 


E , nl , i\ / r i ' r 2 \ ,2m r Za v l /z,m r Zja. _ mr za , . , 

(21 + l)Pi {———)(---r x )*(- T - r 2 ye * 


,2 m r Za 


Xtt 


1=0 


r\r 2 


X 


X 


E 

s=0 


L 2J+1^2m^ ri )Lf +1 (^f^r 2 )s\ 
(s + l + 1 — A)(s + 21 + 1)! 


(4.8) 


Then G' nl (r 1 ,r 2 ) is just the Coulomb Green function evaluated at A = n + 5 A, and taking 
out the pole. In the case where the potentials that appear in Eq. (14.511 are only functions 
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of the modulus of r (i.e. they are rotational invariant), the sum over l reduces to the single 
term that matches the angular momentum l of the bound state. 

Obviously a similar discussion applies to higher order corrections from perturbation the¬ 
ory, and also similar expressions follow for the Lamb shift. 

We will now study each relevant contributing term separately, both in the 1 /m M and 
in the a expansions. We will write explicitly the Z p , Z, pi Z dependence except for the 
dependence on Z that appears in the combination m r Za/m e in the numerical integrals we 
perform. Therefore, such numerical values will change for different muonic atoms. 


4.1 Corrections from the static potential: 

y(°.2) 

4.1.1 One-loop Vacuum Polarization: 5E L VP ~ 0{m p or) 

The Lamb shift in muonic hydrogen, unlike in hydrogen, receives its most important con¬ 
tribution from the electron vacuum polarization. This is due to the fact that the typical 
atomic momentum of the muonic hydrogen is m^a, which is of the order of the electron 
mass: m^a ~ 1.5m e . This effect comes from the modification of the photon propagator, as 
we have already seen in the previous chapter (see Fig. [□). In order to compute it, we must 
take the first order in a of the expansion of II (—k 2 ). 

The integral in r and x can be done analytically. The result reads (see for instance m 

= J d 3 rVyp 2 \r)(p2 P - P2S> = (4.9) 


® / ry \2 

= —( Za) m r 

71 


8it/3 3 1 - 26/S 2 + 352/S 4 - 768/S 6 


4/S 4 (15 - 80/S 2 + 128/S 4 
3 (1 — 4 /S 2 ) 5 / 2 


In 


18(1 - 4/S 2 ) 2 

->/i 


4/S 2 


2/S 


= m r a 3 Z 2 0.005555 


205.00737 meV, 


where 

777/ 

/S = e . = 0.7373836 . (4.10) 

(Zam r ) 

Eq. (14.9(i gives the first entry in Table [1] 

For the case m e <C Zam r the computation can be checked with the result of heavy 
quarkonium. We have checked it. We also observe that m e Zam r is a bad approximation 
to this quantity, so we will not consider it further but only for checking. Actually, neither 
the /3 1 nor the /3 ^ 1 give a good approximation to Eq. (14.9(1 . 


V (0 ’ 3) 

4.1.2 Two-loop Vacuum Polarization: (9(m r .o: 4 ) 

We now compute the 0(m r a 4 ) contribution associated to the two-loop static potential. We 
obtain the second entry of Table [0 

SE^° p3) = m r a 4 Z 2 0.005599695 = 1.50795 meV. (4.11) 
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It agrees with the result of Pachucki [3D] with 5 significant digits. 

We observe that this contribution is significantly bigger than the one coming from double 
insertions of the leading vacuum polarization potential discussed in the next section. In a 
different context a similar situation has been found in heavy quarkonium physics |3T] . 


4.1.3 Double Vacuum Polarization: 6E^ VP x ' vp ~ 0(m r a A ) 

The second 0(m r a A ) correction is generated by the second order perturbation theory of the 
Vyp 2) potential. Following Eq. (14.5j) and the associated discussion we obtain 


y(0,2) v (0,2) 

fipfj WP x v VP __ 


= ( Za) 2 m r ^ 


(v>„iikp 2) 



( 0 , 2 ) 


h)' VYP 


\^nl) 




dp\dp\ 


u (pf) u ( pi ) 
p\ pI 


( [1 + (3{pi + p2)\ 3 f3 2 Q fd 2 {p\ + pi + 2(/3pip 2 ) 2 ] , 

V 12(1 + fa)*{ 1 + fa)* (1 + P Pl )\l + M) 4 


f (1 + /3pi)(l + /3p 2 ) 
V 1 T /3( P i T P 2 ) 


(4.12) 


where 


Q 


12/3 8 PiP2+12/3 8 PiP2 + 36/3 7 p?P2 + p\p A + 12/3 6 p?/>2 + S4/3 6 p\p\ + 74/3 6 p^2 
33/3 5 PiP 2 + 39,3 5 piP2 — 62/3 5 piP2 + 9 f3 A p\ + lllj3 A p\p 2 - 33j3 A p\p\ - l42j3 A p\p\ 

24 pp\ + 99/3 3 pfp 2 - 189/3 3 p 3 p 2 + 18 f p\ - 3p 2 p\p 2 - 75p 2 P 2 lP 2 2 - 24/3p 2 p 2 - 3 p\ 
(pi^p 2 ). (4.13) 


This expression corrects several mistakes in Eq. (11) of [32] (which however gets the correct 
numerical result) and yields the 3rd entry of Table [0 

,,(0,3) vl/(°’ 3 ) /o\ 2 

5El vp VP = m r a 2 Z 2 [-) 0.0055304 = 0.150897 meV. (4.14) 

This numerical value agrees with [3D] within the significant digits given in this reference. 


4.1.4 Static potential (vacuum polarization): 5El ~ 0(m r a 5 ) 

The first four entries in Table Q] are the contributions to the Lamb shift associated to the 
electron vacuum polarization corrections to the static potential V®. Specially difficult is 
the 4th entry, as it corresponds to the three-loop static potential and to the third order 
computation in perturbation theory. It was computed (numerically) in [25] (see also [33] for 
a small correction). It can be split into the following contributions: 

SE^ = m r a 5 Z 2 0.002694 = 5.295 x 10“ 3 meV. (4.15) 


The contribution from 2nd order perturbation theory yields [26] (this result includes all 
permutations): 


SE 


v (°,2) (O.S) 

*vp xi Vp 


+ SE 


V ( 0 ’ 3 ) xV (0 ’ 2) 
v \p x V VP 


= m r a 5 Z 2 0.00109562 


2.153 x 10“ 3 meV. 


(4.16) 
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And the contribution from 3rd order perturbation theory reads [331125]: 


5E 


v (0,2) yV (0,2) v (0,2) 
v yp A v VP A v VP 


= m r a 5 Z 2 0.0000377 = 0.0741 x 10"AneV. 


- 3 . 


(4.17) 


The sum of the above three terms gives the final contribution: 

^Xtic VP = 0.11868m r Z 2 a 2 0) 3 = 0.00752 rneV. (4.18) 

which is the 4th entry of Table HJ The computation has been done independently for a 
time-like (g 2 > 0) and a space-like ( q 2 < 0) momentum of the photon. This last one involves 
the integration of the whole vacuum polarization function II(g 2 ) to the desired order, and 
the other involves just its imaginary part evaluated at g 2 = tm 2 e . 


4.1.5 Static potential (light-by-light): 8E^ ~ 0(m r a 5 ) 

The 5th entry of Table [T] corresponds to the contribution associated to the light-by-light 
corrections to the static potential lA°q i.e. to lAbL (see Eq. (I3.20p ). It was obtained in [27]. 
where a very long explanation was made to argue that the light-by-light contributions could 
be computed in the static approximation. This is evident in the EFT, as they correspond to 
a correction to the static potential, as already stated in Ref. [8]. 

The result for this contribution, given in [27], is 

SE^ = —m r a 5 Z 2 10” 3 [(Z 2 + Z 2 )0.5185 - Z p Z M 0.5852] 

= -m r a 5 0.000452 = -0.00089(2) meV. (4.19) 


4.2 Corrections from the l/m M potentials without vacuum polar¬ 
ization 

We jump directly into the W 2 ) potential, since we already discussed that the V ^ potential 
produces corrections of, utmost, 0(m r a 6 ) and are then beyond the accuracy of our interest. 

We now compute the corrections to the energy and Lamb shift associated to the potentials 
in Eqs. (I3.25|) and (I3.26j) to 0(m r a 5 ). In other words, we compute the relativistic corrections 
that are not associated to the vacuum polarization. 

4.2.1 Relativistic corrections: 5El ~ 0(m r a A ) 

Eq. 113.25(1 is the EFT generalization of the Breit potential. Note that it is in this potential 

m 2 

where the hadronic corrections arise at 0{m r or-^) (we will consider them in more detail 
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later). The energy shift associated to this potential rea> 
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(4.20) 


where 


c 7 .1 2 


Z + l 
’2 z - 1 


3 = 1- 1, 

J = z, 
i — z +1 , 


ha = 


dji-i 


27+3 

— (z + 1 ) j = i - 1 , 
-1 j = z, 

l j = l + 1 , 

— (z + 1 ) ii = z — |, 

3\ = 1 + l ■ 


(4.21) 


(4.22) 


(4.23) 


The energy has been expressed in terms of the total angular momentum J = L + S (where 
S = + S p ) and in terms of the angular momentum of the muon J M =L + S M . The basis 

is taken in terms of the lightest particle, since it is the most convenient one to express the 
energy shift. This is so since the lightest particle gives rise to larger effects in the terms 
which involve the ratio of the masses, and this comes out more clearly when using this basis. 

y(2) 

SE encodes all the 0(m r a A ) corrections to the spectrum due to the 1/m^. It also 
includes higher order effects through the O(a) terms in the NRQCD Wilson coefficients. If 

'In the last line of this equation we have still included the contribution associated to the tau vacuum 
polarization. As its numerical effect is very small we will neglect it in the following. 
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we set them to zero, we obtain the non-trivial leading-order contribution: 
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(4.24) 


We shall also take into account the correction to this order in a coining from the per¬ 
turbative expansion of the relativistic kinetic term, i.e. from Eq. (13.321) . which leads to the 
energy shift: 


SE. 


y(3,o) 

nl 


m 3 Z A a A 
2 m 2 

A 6 



4n 4 n 3 (2l + 1) 


(4.25) 


Summing up the contributions of Eqs. (14.241) , (14.251) we get for the transition of the 
2 S ' 1 / 2 -> 2 P + 2 states: 


8E 


V(2.1) 


5E 


y(3.o) 


m 3 a A Z A m 3 a A Z 3 Z„ 


48 m 2 


16m 2 


4 2 2 

^ r p m p 


= 0.05747- 5.1975 


fm 2 


meV. 


(4.26) 


The first term agrees both analytically and numerically with the one obtained in |30j . 

2 

We shall remark that it has an extra -j+r suppression factor, which was to be expected since 

this correction does not contribute for the case of the hydrogen (in the infinite proton mass 
limit). The 2nd term is the leading contribution associated to the proton radius. Both 
contributions appear as the 6 th and 11th entries in Table Q] 


4.2.2 Relativistic corrections: 5El ~ 0(m r a 5 ) 

We now compute the 0(m r a 5 ) contributions to the spectrum with no electron vacuum 
polarization. As we have already mentioned, this is a well defined quantity, as it amounts 
to the corresponding evaluation of the muonium (fie) spectrum (if we turn off the hadronic 
effects). Taking the 0(m r a 5 ) corrections generated from Eq. (14.201) (typically generated by 
the 0(a) corrections of the NRQCD Wilson coefficients) plus the energy shift produced by 
the expectation value of Eq. (I3.26p . we obtain (note that this computation has been done in 
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the MS scheme) 
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(4.27) 


where /c(n) = ^ is the n-th harmonic number. Note that in this expression the hadronic 

corrections that scale as a 2 : c| ad , C 4 ad and ad 2 ad are also included, as they also produce an 
m r a 5 energy shift. 


4.3 Ultrasoft effects: bEp ~ 0{rn r a 5 ) 



Figure 8: Correction due to ultrasoft photons. 

The interaction of the bound state with ultrasoft photons (symbolically pictured in Fig. [HJ) 
produces an energy shift of 0(m r a 5 ). It has been computed in the MS in Refs. [221 23] for 
the case of hydrogen and positronium respectively. The application to muonic hydrogen is 
trivial, as we only have to rescale for the reduced mass. On top of that we introduce the 
changes for the case of particles with general charges Z M , Z p . Finally, the energy shift reads 
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(in the MS scheme) 
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(4.28) 


where |0 n (O)| 2 = 4 ( mr ^ a ) 3 . In R(n,l) are the Bethe logarithms and are implicitly defined 
by the equality with Eq. (j4.28j) . For their numerical values for the 2 S and 2 P states we have 
used the values quoted in [3D] . 

We observe that SE^ is factorization scale dependent. Such dependence cancels with 
the scale dependence of Eq. (14.271b The sum of both contributions gives all the 0(m r a 5 ) 
corrections to the spectrum that are not associated to the electron vacuum polarization: 


r 77 iO(a 5 ),no—VP 
5E nln „ 


v (2,2) 

SE r&r 


+ 5E^ S 


J nl 


(4.29) 


and is independent of the factorization scale. It can also be split into the different hadronic 
contributions, associated to the fact that the proton is not point-like, and the 0(a 5 ) contri¬ 
bution to the spectrum of two point-like particles (relevant for muonium) in the following 
way: 


rpO(o 5 ),no-VP _ ,pO(a 5 ),no—VP 
^nljj^ point—like 


Ap 0(« 5 ), no—VP 

+ 


5E 


0(a 5 ), no—VP 
nljj^,c^ ad 


+ SE°^TJ yp . 


(4.30) 


Similar equations follow for the Lamb shift energy splitting: 5E^ a f no_VP , although in 
this last case the contribution proportional to c\ &A vanishes, since the spin-dependent term 
does not contribute to the average energy over polarizations. 

The above computation keeps the complete proton and muon mass dependence. It is 
interesting to see the convergence of the m p /m p expansion. We do so for <5-^^^int—nke P , 
which has a non-trivial dependence on this ratio. We obtain 


0(m p a 5 ) : 

rnO(a 5 ),no—VP _ 
point—like 

—0.900 meV 

o(V|) : 

A pO(tt S ),no-VP _ 
point—like 

0.226 meV 

0(m p a 5 E |) : 

Ap O(a 5 ),no-VP _ 
point—like 

—0.054 meV 

0(m p a^) : 

AF 0(« 5 ),no-VP _ 
point—like 

0.010 meV, 


which, summing up to all orders, leads to the following 0(a 5 ) energy contribution to the 
Lamb shift 

7 = -0.71896 meV, (4.31) 
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which corresponds to the 7th entry of Table [Lj This result is very similar to the one computed 
by Pachucki [301, where these effects sum up to U(2 Pi/ 2 — 2 S 1 / 2 ) = —0.663 — 0.045 — 0.010 = 
—0.718meV at C>(m M a 5 ). 

We now consider the hadronic corrections. The energy shift associated to the hadronic 
vacuum polarization reads 

^ = 0.0111(2) meV, (4.32) 

which corresponds to the 14th entry of Table [Q 

The energy shift associated to C 3 ad is usually named the two-photon exchange contribu¬ 
tion. Using the lower value in Eq. 02.221) we obtained 

\ = 0.0344(125) meV, (4.33) 

which corresponds to the 15th entry of Table [TJ 
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4.4 1/mj) electron vacuum polarization corrections: 5 El ~ 0{m r a h ) 

We now compute the energy shifts, with 0(m r a 5 ) precision, associated to the electron vac¬ 
uum polarization. They are produced by 2nd order non-relativistic quantum mechanics 
perturbation theory of Uyp 2 ^ ~ a 2 /r, together with the V^ 2,1 ) ~ a/m 2 and W 3,0 -’ ~ 1/m 3 
potentials, as well as by the correction due to the Uyp 2 ^ ~ a 2 /m 2 potential. This sum con¬ 
stitutes a well defined set, as it can be parametrically distinguished from other contributions 
(formally through the number of light fermions). The energy shift then reads 
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For the Lamb shift corrections we obtain the explicit expressions 




(4.34) 
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= ^0.046473 — 0.016953^2^ meV + C>(a 6 ). 

For this last result we have used Eq. (14.8ft . Summing up both contributions, 
(I4.36p . we obtain 

y(2.2) y(2,l) x y(0,2) y(3,0) x y(0,2) 

5El VP +5 El VP +6 El VP 

rp 2 / j,2 

= m r a 5 0.0095460 - m r a 5 0.01433-^ = 0.018759 - 0.0282^7 


(i + q py 


(4.36) 

Eqs. (I4.35P and 


As we have already stated, this sum constitutes a well defined set, as it can be parametrically 
distinguished from other contributions (formally through the number of light fermions). This 
is also so for each individual term in the last equality in Eq. (I4.37p . The first term corresponds 
to assuming the proton to be point-like (switching off the proton radius contribution) and 
gives the viii) entry in Table [[] This contribution was first computed in pjj and later 
corrected in [3T, [35] - Nevertheless, a different number has been obtained in two recent 
analyses |29) 36]. We confirm this last number, which is the one we quote in Table [0 

The term proportional to the proton radius in Eq. (I4.37p emanates from the coefficient 
U. It corresponds to the xi) entry of the table and it is in agreement with the result found 
in [30]. 


4.5 0(m r a 6 x In) effects 

The first 8 entries in Table[Qgive the complete 0(m r a b ) result for a point-like proton. In this 
result we have kept the exact mass dependence. The 0{m r a 6 ) contribution is dominated by 
the logarithmic enhanced terms. Here, we compute the leading ones. We assign a general 
counting of m r /m p < In a ~ In(m e /m M ). Therefore, we only compute those contributions 
at leading order in the m r /m p expansion, i.e. those where the proton is infinitely massive. 
If we switch off electron vacuum polarization effects (i.e. we switch off the interaction with 
the electron) the system corresponds to the standard hydrogen situation, which has no 
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0(m r a 6 In a) effects. Actually, this is also true if we consider the case of muonium (with 
finite recoil effects), which again has no 0(m r a 6 In a) effects. The reason is that the sum of 
all possible contributions vanishes for the case of the lamb shift, since the effective energy 

shift is [370 

,5 


$ E n i s 


m; 


3 2 
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6 In — 

a 


_ 7 17 ) 


n° 


(4.38) 


which vanishes for the Lamb shift. Therefore, we can actually claim that all the 0(m r a e In a) 
logarithms are generated by the electron vacuum polarization (for a point-like proton). Note 
that this would also be true if we consider proton-recoil corrections. In any case, as we have 
already mentioned, here we only consider the infinite proton mass limit. In this limit, for a 
point-like proton, only two contributions are produced (both of them generated by electron 
vacuum polarization effects), listed in the the ix) and x) entries of Tabled! which we now 
discuss. 

The 9th entry is due to the logarithmic-enhanced 0(a 2 ) corrections to the c'^ Wilson 
coefficient (see Eq. (12.10)) ) that appear in the tree-level potential (see the c^-dependent term 
of Eq. (13.251) 1. It produces an ct 3 /m^xlogarithm-potential, the expectation value of which 
gives the following energy shift to the spectrum 


5E nl = 


S 4 (m) 

m*a 


2 rrif. 


n° 




0(a Q In) 


(4.39) 


and to the Lamb shift 


5E l = -m r a e 0.08885 = -0.0012741 rneV, (4.40) 


which is the number that we quote in the 9th entry of Table [T] 

The 10th entry in Tabled] is generated in the same way as the 8th entry but multiplied 
by the (logarithmic enhanced) O(a) term of Cp\u) (see Eqs. (12.10)) and (14.35)) ): 


f°° j 2 u ( ( 1 2 ) c d ) f(l + 2 /^)(l + 2/^(l + /3g)) 

[mra) 16 7T J 4 q q 2 ml\ (l + /3g) 4 

4(3Y + 12/3 3 (? 3 + 4 (3 2 q 2 + ll/3q + 3 4 + X ) ln (^Ti) 

(Pq + 1) 5 (/3g + l) 4 


(4.41) 

G(a Q ln) 


The v dependence gets regulated by the ultrasoft scale, which we set to v = m^a 2 , producing 
the number 


v (2,2) 

5El VP 


= -m r a 6 0.31644 = -0.004538 rneV, 


(4.42) 


0(a Q ln) 

which we quote in the 10th entry in Table [Q 


s For simplicity we set Z p = = 1 in this section. 
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Both computations were considered before in Ref. |30| . We agree with them for the 
significant digits given in that reference, ft is also interesting to see that both contributions 
can be understood from a renormalization group analysis in some appropriate limit [38] . 
This analysis also gives us information on the logarithmic structure of the recoil, m r /m p , 
corrections. At this order extra logarithmic terms appear. Nevertheless, they are at most 
linear: (9(m M a 6 — In a), i.e. there are no 0(m r a e — In 2 a) terms, contrary to the claim in 
Ref. [29]. 

For a point-like proton this computation would finish our analysis. The error would be 
due to uncomputed contributions of 0(m r a 6 ) and 0(m r a 6 ^ In a). In Refs. ISjES] several 
terms of this order were computed. We use these analyses to estimate the error. Specially 
useful to us are the (a) and (d) entries in Table IV of the last reference. They are related with 
the large logarithmic contributions discussed above but also include some finite pieces. We 
take the difference with the pure logarithmic terms for the generic 0{m r a 6 ) error. Taking 
instead 1/2 of the sum of the 9th and 10th entries yields a similar error: ~ 3 /ieV. This is the 
error we quote in the first term of Eq. (15.11) . which encodes all the QED-like contributions 
assuming the proton to be point-like. 

We now consider the 0(m r a d In a) correction associated to the proton radius. It scales 
like 0(m r a G In a x ml r 2 ) and has been computed in [40]. Such effect would be generated by 
the 2nd order perturbation theory of the delta potential (note that a similar effect would also 
exist in the analogous hydrogen computation). The infrared behaviour of this computation 
would be regulated by the inverse Bohr radius generated by the bound state dynamics, 
~ m p a. The ultraviolet behaviour gets regulated by energy scales of order ~ m T . This 
produces the large logarithm: In {{m p a)/m p ) = In cc. The explicit correction reads 


SE l 


2tt 

T 


raj; or 

2 % 


r 2 v c? In a 


-0.0014 


fm 2 ’ 


(4.43) 


and it is listed in the 13th entry of Table [0 We use 1/2 of this result for the error of the r 2 
coefficient in Eq. (II.2|) and Eq. (15.1(1 . 

A complete analysis of the 0(m r a 6 ) effects from an EFT perspective will be discussed 
elsewhere. 


5 Summary of results and conclusions 


All the contributions to the Lamb shift considered in this paper are listed in Table [T] Their 
sum produces the following theoretical prediction for the Lamb shift 


^ithiswork 


206.0243(30) 


5.2270(7)-^r + 0.0455(125) 
fm 2 


meV. 


(5.1) 


The first ten terms in Table Q] are those associated to a pure QED-like computation assum¬ 
ing the proton to be point-like. Their sum is the first term in Eq. (15. ip . and its error is the 
estimate of the 0(m r a 6 ) effects. The second term in Eq. (15.ip encodes all the corrections 
proportional to the proton radius x)-xiii) entries in Table [0. The error of the coefficient 
of the term proportional to r 2 is the estimated size of the (9(m r a 6 (m 2 r 2 )) terms. The last 















term encodes the r 2 -independent hadronic effects. The error is the assigned uncertainty due 

m 3 

to unknown terms of 0(m r a 5 -^). We emphasize that a partial incorporation of subleading 
corrections in a to the above expression will not improve the precision of the result (unless 
there are arguments to argue that such contributions are the dominant ones), as the uncer¬ 
tainty is still dominated by unknown parametric terms of order m r a & . For an account of 
some of these corrections see [41] . 

In order to obtain our primary result Eq. (II.2|) . the first and last term of Eq. (15. ip has 
been added and the error combined in quadrature. With this, together with the experi¬ 
mental result in Eq. «LU), we obtained in Ref. [6] the value for the proton radius quoted in 
Eq. (II.3j) , where the theoretical and experimental errors have been combined in quadrature. 
Nevertheless, the latter is completely subdominant with respect to the total error, which is 
fully dominated by the hadronic effects. In this respect it is also convenient to present our 
result in the following way 


a 


m 2 v fm 2 


A E l = 206.0243 meV 
1 m 3 a 3 

71 8 

1 m 3 a 3 

71 8 

+0(m r a 6 ) 


[47.3525 + 35.1491a + 47.3525a 2 ln(l/a)] 


-4 [co ad + 16vra4 ad l 

my 3 J 


(5.2) 


Note that since Cg ad ~ a 2 and ad if" 1 ~ a 2 , the third line of the previous equation encodes all 
the hadronic effects that are not related to the proton radius of order a 5 . This presentation 
of the result where r p and c| ad are kept explicit could be important for the future. In the 
long term (once the origin of the proton radius puzzle is clarified) the natural place where 
to get the proton radius is from the hydrogen Lamb shift and Cg ad (once the radius has been 
obtained) from the muonic hydrogen, since Cg ad is suppressed by an extra factor of the lepton 
mass. In this scenario a complete evaluation of the 0(m r a e ) term may improve the precision 
of an eventual experimental determination of c| ad . Note that in this discussion we assume 
that we can determine d\ &d from alternative methods, like dispersion relations. 

Finally, we profit this computation to give in the Appendix the exact a 5 expression for the 
muonium spectrum, keeping the complete mass dependence, which can be easily deduced by 
changing m p m p and m p —$■ m e , and setting the hadronic coefficients, dj, and the electron 
vacuum polarization effects to zero. 
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i) 

0(m r a 3 ) 

y(°) 

V YP 

Eq. (SD 

205. 00737 

ii) 

0(m r a A ) 

T/W 

V VP 

Eq. (14.111) 

1. 50795 

hi) 

0(m r a 4 ) 

T/W 

V VP 

Eq. (14.141) 

0. 15090 

iv) 

0(m r a 5 ) 

V VP 

Eq. (14.181) 

0. 00752 

v) 

0(m r a 5 ) 

v LbL 

Eq. (14.191) 

-0. 00089(2) 

Vi) 

°(m r a 4 x ^|) 

y(2,i) _)_ y(3,o) 

Eq. (14.261) 

0. 05747 

vii) 

0(m r a 5 ) 

14 o-vp + ultrasoft 

Eq. (14.311) 

-0. 71896 

viii) 

0(m r a 5 ) 

V™ + E (2,1) x l4°p 2) + • • • 

Eq. (I4.37P 

0. 01876 

ix) 

<D(m r a 6 x ln(^)) 


Eq. (14.401) 

-0. 00127 

x) 

0(m r a 6 x In a) 

_ Ui) 

K VP : C D 

Eq. (14.421) 

-0. 00454 

xi) 

0(m r a A x mfr 2 ) 

y(2,l) (P) 2 

v i c D i ' p 

Eq. (14.261) 

-5. 19745^ 

xii) 

0(m r a 5 x mfr 2 ) 

T/( 2 ’ 2 ) r (P). r 2 

v VP ' >>'p 

Eq. (I4.37P 

-0. 02815^ 

xiii) 

0(m r a 6 In a x vnfr^) 


Eq. (I4.43P 

-0. 00136^ 

xiv) 

0(m r a 5 x g) 


Eq. (14.321) 

0. 0111(2) 

xv) 

0{m r a 5 x 

\ ' m~ m 7T / 

y(2). c ha d 

Eq. (|4.33p 

0. 0344(125) 


Table 1: The different contributions to the Lamb shift in muonic hydrogen in meV units. 

A Muonium spectrum 

We profit from the results obtained in this work to give the spectrum for the muonium bound 
state (/re) for general quantum numbers at 0(m r a b ). We first exchange the proton by the 
muon and the muon by the electron. Then, the main difference with muonic hydrogen is 
the lack of hadronic contributions, as well as the fact that all electron vacuum polarization 
effects can be eliminated, in particular this implies that the static potential becomes trivial. 
Thus, we are only left with the relativistic corrections to the potential which come from Eqs. 
(I3.2ip and (I3.22[) plus the energy coming from the kinetic term and the ultrasoft effect. The 
ultrasoft correction to the energy only depends on the reduced mass, and so it will be the 
same as the one for the muonic hydrogen in Eq. (I4.28p . Altogether, for a given energy level 
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(A.l) 


where Cjj, hjj and dj e j have been defined in Eqs. (14.2ID - (14.231) . and the first and second 
parenthesis in the right hand side of the first equality contain the 0(m r a 4 ) and 0(m r a 5 ) 
contributions respectively. Note that the exact mass dependence has been kept in this 
expression to order a 5 . 

The expressions for the potential of muonium can also be found in Ref. 02. One could 
be worried that the potential is different to the one we use. The reason for this difference is 
that they obtain the potential by matching on-shell S-matrix elements (and by a change in 
the renormalization scheme of the ultrasoft computation), still their potential is equivalent 
to ours through field redefinitions, and yields the same physical results. In particular, for 
spin-independent states the result for the energy shift can already be found in Eqs. (2.12) 
and (2.13) of that reference. 
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